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Abstract: The quasar 3C 279 is one of the best-studied flat spectrum radio quasars. It is located at a 
comparatively large redshift of z = 0.536: E > 100 GeV observations of such distant sources were until 
recently impossible both due to the expected steep energy spectrum and the expected attenuation of the 
y-rays by the extragalactic background light. Here we present results on the observation of 3C 279 with 
the MAGIC telescope in early 2006. We report the detection of a significant very high energy y-ray signal 
in the MAGIC energy range on the observation night of 2006 February 23. 
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Introduction 

3C 279 (RA=12h56''ll.P, Dec=-5°47'"22') was 
the first blazar discovered in y-rays with the Comp- 
ton Gamma-Ray Observatory [1]. It is an excep- 
tionally bright and variable source at various en- 
ergy bands. In the optical, typical variations of 
2™ have been observed, while during flares, up to 
8™ were reached [2]. Variations were observed on 
day-scale [3], are mostly coherent at all frequen- 
cies, and the spectral energy distribution (SED) 
hardens when the source brightens [4]. Typical 
flux variations by factors of ^ 20 in the GeV and 
~ 5 - 10 in the IR-to-UV have been found [5]. The 
strongest variabiUty occurs on timescales of a few 
weeks to ~ 6 months. 

Blazars are thought to be supermassive black 
holes in the centers of galaxies accreting matter 
They possess strongly collimated, ultra-relativistic 
plasma outflows (jets), aligned closely to the ob- 
server's line of sight. Thus, their SEDs are al- 
most entirely dominated by the jet emission. In 
v/v representations, they are characterized by two 
distinct nonthermal components. The lower en- 



ergy bump is commonly ascribed to synchrotron 
radiation emitted by relativistic electrons. There 
is less agreement about the origin of the high- 
energy bump. In leptonic acceleration models, 
inverse Compton (IC) scattering of synchrotron 
(self-synchrotron Compton, SSC, models) or am- 
bient photons (external-inverse Compton, EIC, 
models) on high energy electrons explain the MeV- 
to-TeV radiation. While BL Lac objects can often 
successfully be described by SSC models, the more 
luminous Flat Spectrum Radio Quasars (FSRQ) 
usually are modeled requiring external compo- 
nents, particularly in the £ > 100 keV regime [6], 
although the existence of EIC components is ques- 
tioned e.g. by [7, 8,9]. 

Since its discovery in y-rays, 3C 279 was exten- 
sively studied, in particular also during various 
multi-wavelength campaigns [3, 4, 6, 10, 11, 12]. 
The broadband SED extends from radio frequen- 
cies to the y regime and is, with a comparatively 
low synchrotron peak, fairly typical of FSRQs. 
Correlations between the two peaks, which are ex- 
pected in leptonic models, were partially observed 
[6], although no consistent patterns were found. 



Discovery of VHE Gamma-Rays from 3C 279 



The SED of this source is rather complicated, and 
in spite of detailed observations, still poorly under- 
stood. Observations of the high-energy part of the 
spectrum are complicated by the fact that terres- 
trial and satellite-borne instruments together can, 
at present, not fully cover the frequency range of 
the high-energy bump. Additionally, the very high 
energy (VHE, defined hy E> 100 GeV) spectrum 
is in part suppressed from interactions with the ex- 
tragalactic background light (EBL). This modifica- 
tion is, of course, strongly dependent on the source 
distance. Thus, distant VHE 7-ray sources repre- 
sent an excellent tool for determining the 0.3 to 
30 nm EBL (e.g., [13]), which at redshift z = 
was observed by various satellite experiments, al- 
though direct measurements suffer from huge fore- 
ground contaminations by light contributions from 
the solar system and our galaxy. As of now, sev- 
eral different EBL models have been proposed 
[13, 14, 15]. A precise measurement of the en- 
ergy spectrum of 3C 279 is crucial for two reasons: 
With a state-of-the-art EBL model, emission mod- 
els for 3C 279 can be tested in detail. By using con- 
servative arguments on particle acceleration mech- 
anisms, the y-ray emission of 3C 279 also permits 
to formulate stringent constraints on the EBL level. 

All extragalactic VHE y-ray sources detected so 
far are of the BL Lac type. These objects are the 
low-luminosity counterparts of the FSRQ class of 
AGNs [16], with their synchrotron peaks shifted 
to higher energies [17]. BL Lac objects have 
been detected aplenty recently in the VHE range: 
17 BL Lacs have been found so far,^ reaching 
to redshifts of z = 0.212 [18], plus M87, as- 
sumed to be a misaligned BL Lac. The sam- 
ple includes also blazars as PG1553-H113 [19] 
with its extremely soft energy spectrum and re- 
cently discovered BL Lacertae itself [20], which 
is the first "low-peaked BL Lac object" (Vsynchr < 
10'^ Hz). Both observations were largely made 
possible from the exceptionally low energy thresh- 
old of the MAGIC telescope. 

In January 2006, 3C 279 was found in a high op- 
tical state, brightening to 14™5 in the R-h&nd} A 
WEBT campaign in early 2006 [12, 21] and multi- 
wavelength campaigns during both periods [11] 
were performed in order to get further information 
on the temporal and spectral properties of 3C 279. 



Up to recently, a VHE y-ray detection of 3C 279 
was prevented by its high redshift of z = 0.536. 
The resulting cutoff due to EBL attenuation is ex- 
pected at around E « 200 GeV. The MAGIC (Ma- 
jor Atmospheric Gamma-ray Imaging Cerenkov) 
telescope [22] is currently the largest single-dish 
(17 m diameter) Imaging Atmospheric Cerenkov 
Telescope, located on the Canary Island of La 
Palma. With its low trigger threshold (50 — 60 GeV 
at low zenith angles) is best suited for observing 
the lower part of the VHE range of distant AGN 
spectra. At the geographical latitude of MAGIC 
(28°45' N), 3C 279 can be observed under medium 
zenith angles (above 34°), with an accordingly in- 
creased observation threshold. 

Observations and Analysis 

3C 279 was observed from late January to April 
2006 and in January 2007. Simultaneous opti- 
cal /J-band observations were carried out using the 
1.03 m telescope at the Tuorla Observatory, Fin- 
land, and the 35 cm KVA telescope on La Palma. 
Here we report results from the 2006 observations. 

The VHE y-ray observations were performed in 
the ON-OFF mode. ON data were collected while 
pointing directly to the source, while OFF data, 
necessary for the background estimation, were 
recorded by pointing to a nearby region of the sky. 
The OFF region, in which no y-ray source is ex- 
pected, was chosen to have a comparable zenith 
angle distribution and night sky conditions. From 
2006 January 29 to 2006 April 4, 3C 279 was ob- 
served for 14.9 hours. In addition, 3.9 hours of 
OFF data were recorded. Data taken during non- 
optimal weather conditions or affected by hard- 
ware problems were excluded from the analysis. 
This concerned 5.2 hours worth of ON data. The 
remaining events were calibrated [23] and ana- 
lyzed using the standard MAGIC analysis chain 
[24, 25 1. Briefly, the analysis proceeds by reduc- 
ing the image to a single light cluster by removing 
noise, calculating image parameters [26], and us- 
ing a multi- variant method to discriminate y-like 

1 . See http://www.mppmu.mpg.de/~rwagner/sources/ 
for an up-to-date list 

2. Long-term optical monitoring data are available at 
http://users.utu.fl/kani/lm/ 
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events from the dominant hadronic background. 
For the latter step, we use the Random Forest (RF) 
method [27], which combines the image parame- 
ters into a single test statistic, called hadronness, 
based on training samples of both Monte Carlo 
generated 7 [28, 29] and real background events. 
The excess events were identified using the classi- 
cal ALPHA approach by subtracting suitably nor- 
malized OFF data from the ON events. The sig- 
nificance of any excess was calculated according 
to eq. 17 in [30]. The energy estimation for each 
event was performed using a RF method, too, lead- 
ing to a reasonably constant energy resolution of 
«23% above 150 GeV. 

Results 

The data were subjected to three independent anal- 
ysis chains, which obtain compatible results w.r.t. 
each other A standard analysis with an energy 
threshold of w 200 GeV was used to infer the E > 
200 GeV light curve given in Figure 1 . Also shown 
is a 7? -band optical light curve for ten MAGIC ob- 
servation nights from 2006 January 31 to 2006 
March 3 1 . The typical observation time per night is 
around one hour. While in most of the nights 7-ray 
fluxes compatible with zero were observed, dur- 
ing the 2006 February 22 observations a marginal 
signal was seen. In the night of 2006 Febru- 
ary 23 we found a clear 7-ray signal with an in- 
tegrated photon flux F{E > 200 GeV) = (3.5 ± 
0.8) X 10^"cm^^s^^. The source was observed 
for 62 minutes (MJD 53789.1633-53789.2064) at 
zenith angles between 35° and 38°. The stable 
event rate during this observation allows classify- 
ing it as "dark night" rate, although minimal moon- 
shine was present. 

A low-energy analysis, exploiting the timing prop- 
erties of the air shower images with a lower anal- 
ysis threshold, « 110 GeV, was used to calculate 
the significance of the found excess. Threshold is 
defined as the maximum of the energy distribution 
of the accepted events, viz. showers with ener- 
gies down to w 80 GeV are included in the sam- 
ple. Both analyses show consistent results from 
E > 200 GeV on. After subjecting the ON data to 
an ALPHA cut inferred from training on an (inde- 
pendent) Crab nebula data sample recorded at sim- 
ilar zenith angles, and subtracting a properly nor- 
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Figure 1: MAGIC E > 200 GeV 7-ray (top) and 
optical 7?-band (bottom) light curves obtained for 
3C 279 in early 2006. 



malized OFF distribution in the signal region, 624 
and 93 excess events remain between 80 and 220 
GeV and between 220 and 600 GeV, respectively 
(Fig. 2). The data were separated into these two 
independent samples because of the very different 
7/hadron separation powers in the two energy re- 
gions: Adding the highly-enriched, low statistics 
high-energy 7 sample to the large statistics low- 
energy sample would spoil the overall significance 
calculation. The resulting significances are 6.1(7 in 
the low energy region and 5 . 1 (7 in the high energy 
region; the excess is compatible with a point-like 
source and its position is consistent, within statis- 
tical uncertainties, with the 3C 279 position. The 
signal at low energies exceeds the intensity of the 
Crab nebula, while for E > 200 GeV it is «15% 
of the Crab nebula flux. The spectral analysis of 
the data is ongoing. Our 7-ray detection was not 
accompanied by an optical flare or by particularly 
high flux levels or outbursts in X-rays [31, 32]. 

Conclusions 

In observations triggered by a high optical state of 
the flat spectrum radio quasar 3C 279, MAGIC was 
able to detect a highly significant VHE 7-ray sig- 
nal from this source, well-known and much studied 
at lower energies. This discovery is a crucial step 
forward for VHE 7-ray astronomy in various con- 
texts: 
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Figure 2: The 7-ray signal below (left) and above 
(right) 220 GeV on 2006 February 23. Showers 
which are aligned parallel to the telescope axis 
have small ALPHA values. To remove residual 
hadronic background, OFF data were normalized 
to the ON distribution between 20° < ALPHA < 
80° and subtracted from it. The structure visible in 
the left ALPHA plot in the OFF data is statistically 
compatible with a flat distribution. 



• The distance over which astrophysical ob- 
jects can be observed at VHE 7-ray energies 
was enlarged substantially — when the pho- 
tons recorded during our observations left 
3C 279, the age of the Universe was only 8.4 
Gyr. 

• For the first time, VHE 7-rays from a flat 
spectrum radio quasar were detected. These 
objects are very luminous, but are not ex- 
pected to show intense £ > 100 GeV emis- 
sion. With 3C 279 and BL Lacertae, ob- 
jects of all classes comprising the "blazar 
sequence" [17] have now been detected in 
VHE y-rays. 

• A precise measurement of the energy spec- 
trum in the VHE region may allow for strin- 
gent constraints on the EBL density. 
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